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^ ■ Abstract 

cn 

2^ ' In this paper, we study the Euler and Euler-Poisson equations in R^ , with multiple 7-law 

for pressure function: 

, , , vn, 

X' P(p) = e^^p-^ (1) 

. ^ ; 

where all 7^+1 > 7i > 1, is the constants. The analytical line solutions are constructed for 
the systems. It is novel to discover the analytical solutions to handle the systems with mixed 
pressure function. And our solutions can be extended to the systems with the generalized 
multiple damping and pressure function. 
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1 Introduction 

The TV-dimensional Euler equations can be formulated as the follows: 

Pt+V • {pu) =0, 

p [ut + (u • V) u] + '^P{p)= -6p\j^- pF{t), 

(2) 
St + u-VS ^0, 

A$(t,f) =a{N)p, 

with X = (xi, X2, ••■, xn) £ R^ , 

and p = p{t, x), u{t, x) e R^ and 5" — S{t) are the density, the velocity and the entropy respectively. 

And a{N) is a constant related to the unit ball in i?^: a{l) = 2; a(2) = 27r and for iV > 3, 

a(JV) ^ NiN ~ 2)VoliN) ^ JV(JV - 2) ^^^^^ ^ ^^ , (3) 

where Vol{N) is the volume of the unit ball in R^ and F is a Gamma function. 

When S ^ 1, the system can model fluids that are self-gravitating , such as gaseous stars. For 
N = 3, the equations ([2]) are the classical (non-relativistic) descriptions of a galaxy in astrophysics. 
See [I] and [2], for details about the systems. 

When (5 = — 1, the system is the compressible Euler-Poisson equations with repulsive forces. The 
equation ([2|)4 is the Poisson equation through which the potential with repulsive forces is deter- 
mined by the density distribution of the electrons. In this case, the system can be viewed as a 
semiconductor model. See [3] and [10] for detailed analysis of the system. 

When (5 = 0, the potential forces are ignored. The system is called the Euler equations. See [4] 
and [13] for detailed analysis of the system. 

Here P ^ P{p) is the pressure, the 7-law on the pressure for the single gas, i.e. 

Pip) = e>^ (4) 

is a universal hypothesis. The constant 7 ~ cp/cy > 1, where Cp and Cy are the specific heats per 
unit mass under constant pressure and constant volume respectively, is the ratio of the specific 
heats. In particular, the fluid is called isothermal if 7 = 1. More generally, the pressure function 
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of mixed gases, can be expressed by the multiple 7-law ([H], [H], and [TS]): 

m 

P{p) =Pi+P2 + ...Pn = e'J2 P^' ' 



(5) 



where all 7j+i > 7j > 1, is the constant, m in a, positive integer. 

The system with the multiple 7-law can reflect the better approximation of the real situations. 
For example, the fluids in the stars mix with many types of gases in models of gaseous stars in 
astrophysics [T3]. 

Here the time-dependent external force F{t) — {Fi{t), F2{t),..., Fjv(O) G C'° are coupled in the 
systems. 

The system ^ can be rewritten in scalar form. 



dp 
at 



^^' dp 



I V^" op I V-^JV 



dxk ■ 



=0, 



For the Euler equations ©, (6 = 0), with F{t) — 0, in radial symmetry: 



X X 

p{t, x) — pit, r) and u ~ —V{t, r) := —V, 

r r 



ij) 



withr= E^=l 



-^ ^2 



1/2 



there exists a family of solutions for the Euler equations ©, ((5 = 0), without the external force 
(F = 0), for 7 > 1, [9] 

H(r/a(t))l/(T-l) 



p{t,r) = 



a(t)« 

0, 



for 2/(^) > 0; 

for viw)) < 



Vit^^) = W)''' Sit,r)^lnK, 



(8) 



"W = ^7TT+iv(7^TT' «(0) = flo > 0, a(0) = oi. 



y(a;) 



_ (7-l)A 2 I Q,e-i 



where iiT is a positive number, 
for 7 = 1, [H] 



„y{r-/ci(t)) 



p{t,r)^^-^^^,Vit,r)^^^r, S{t,r)^\nK 



"W = ^'«(0) = ao > 0,a(0) = ai. 



(9) 



y(x) - ^-^2 
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where A, a, ao and ai are constants. 

The separation method for analytical solutions, were used to handle other similar systems with 
single 7 functions: 

P{p) = e^p\ (10) 

or without pressure function, in [7], [5], [9], [l2], [16], [H], [H], [E] and [20]. 

It is very natural to extend the results for the system with multiple 7 function ^: 



P{p) = e^Y.' 



ai 



(11) 



j=i 



In this article, we have obtained a class of line solutions to the Euler equations ^ ((5 = 0) and 
Euler-Poisson equations ^ (5 = ±1), with multiple Gamma function, in the following theorems: 



Theorem 1 For the Euler equations with multiple Gamma Junction, [5 — 0) ^, we have the 
family of the solutions, 



N 

p{t,x) = / ( E Ci [xi - ai[t)) ) , 



u{t,x) ^{ai{t), 02(0, •■•, aN{t)), 



S{t, x) ~ In 



N 

g [ T, Ci{xi - ai{t)) 

i—l 



where 



&^{t) ^ F,it) + aC, for I = 1,2,. ..N, 
a(0) = ao, a(0) = ai. 



(l)for 71 > 1 , with 



7i-2 



? + 5(^)E"i/(^P"'+.9WE"i7.7(^r'"V(^)=0 



(12) 



(13) 



(14) 



g{z) > 0, for z e (-00,00), 

where Ci, C2, ■■-, Cn, ^, go O'l^-d ffi are arbitrary constants; and f > is an arbitrary C^ function 
(2) for 71 = 1, with 

i + pC,g{z) + 5(z)fgc. + g{z) Y.%2 f i^V'^' + g{z) ^7=2 7,7 (^)"'"' / (z) = 0, 

g{z) > 0, for z e (-00,00), 

f > is an arbitrary C^ function. 



(15) 
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We notice that the velocity u{t,x) of the solutions ([T^ and (US]), are only time-dependent 
functions: 

u{t,x) = (ai{t),a2{t),...,aN{t)). (16) 



It is different from the conventional velocity for the analytical solutions: 



^ alt) ^ 

a{t) 



(17) 



Moreover, we need some modification to have the corresponding results for the Euler-Poisson 
equations. 

Theorem 2 For the Euler-Poisson equations with multiple Gamma function, (S — ±1), (0), we 
have the family of the solutions, 



^ Cf .. / N \ 

P(t,x) = ^^/ i^C^ {x, - a,{t))\ 

u(t,x) ^ {di{t), 02(0, •■•, a-Nit)), 



(18) 



S{t,x) = In 



N 

g ( E Ci{xi - a^{t)) 



where 



■d,{t) = F,{t) + C\^ + d^{t), for ^ = 1, 2, ...N, 

a(0) — uq, d(0) = oi, 
(l)for 71 > 1, with 

e + 4^9(z) E7=i / izf^-' + ^9(z) EU ^^f (^)''"' / (^) + ^/(^) = 0> 

1=1 ^ 

g{z) > 0, /or z G (-00, 00), 



(19) 



i = l ^ 



(20) 



Af 



where Ci, C2, ..., Cat, ^, go a^c' 5i O'^e arbitrary constants with S C~ > 0; / > is an arbitrary 



=1 



C^ function; and di,d2, ■■■■,d]s[ are arbitrary C° functions; 
(2) for 71 = 1 , with 



a{N) 



"W „f^\ fi^) 



^+^^g(z) + ^^giz) 



+4^9iz) E%2 Hzp-' + 4^9{z) E%2 i^fi^r-' 7 (-) + ^/(-) = 0, (21) 

i = l ^ 

9{z) > 0, for z e (-00,00), 



i = l ' 



f > is an arbitrary C^ function. 
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Remark 3 The mass of the solutions il2\) for the Euler equations, and filf^) for the Euler-Poisson 
equations, in 1-dimensional case, is finite, if 

f{z)dz < oo, (22) 

J —OO 

and 

/oo 
f{z)dz < oo, (23) 

-oo 

respectively. 

2 Line Solutions 

With regard to the continuity equation of mass ©I, we found the following solution structures of 
the below lemmas fit it well: 

Lemma 4 For the mass equation: 

Pt+V ■ (pu) =0, (24) 

there exist solutions, 

p{t, f ) = / ( E Cj {xi - Oiit)) j , u{t, x) =u = (ai(i), a2{t), ..., aN{t)), (25) 

with the form arbitrary f > £ C^ and arbitrary a, (t) G C^ . 
Proof. For the mass equation, we have 

pt + Vp-u + pV -11 (26) 

= ^/ (^a {x, - a,{t))\ + V/ r§Q {x, ~ adt))\ ■ {ai{t), d2(t), ..., djv(t)) (27) 

= / (^a (x. - a,{t))) (-^^CMt)] (28) 

+ / (^a {x^ - a,{t))] (Cidi(i) + C2a2{t) + .... + CNCiNit)) (29) 

= / [ E C, {x, - a^{t)) - E CMt) + S CMt) (30) 



i=l 



(31) 



The proof is completed. 
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Similarly, we have the corresponding lemma for the entropy equation ([2]) 3 : 

Lemma 5 For the entropy equation: 

St + u-\7S = 0, (32) 

there exist solutions, 



N 

S{t,x) = G I Y^ Ci{xi- ai{t))\ , u{t,x) =u = {ai(t), a2{t),..., aAr(i)), (33) 



with the form / > G C^ and ai{t) e C^. 
Proof. For the mass equation, we can obtain: 

St+u-yS (34) 

= ^G Ts C, (x, - a^{t))] + (di(i), a2{t), ..., dw(i)) • VG (^C, (x. - a,{t))\ (35) 
= G (^a {x^ - a,;(i)) j T- J^G,d,(t) j (36) 

+ G LS G, (a;, - a,(i)) j (Gidi(i) + G2d2(t) + .... + GjvdA'(i)) (37) 

= 0. (38) 

The proof is completed. ■ 

The technique of constructing solutions is to deduce the partial differential equations into 
ordinary differential equations only. Based on the above lemmas, it is clear to check our solutions 
for the system. 

Proof of Theorem [TJ Our structure of the solutions (IT^ . fits well for the mass equation ^i 
and the entropy equation ^3, from the above lemmas. For the i-th momentum equation ^2, we 



N 

can define z :— T, Ci {xi — ai{t)) , to have 

i—l 
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N 



M^-E-'i: hi-.'i:.-+.«(') 



dui Y^ 9ui \ d 

f-J *" dxk j dxi 

d 



oxj ^-^ 



j=i 



- P«(t) + ( ^9 {z)\ aJ2f i^V + 9{z)^ I Y. f (^)^' I + pp^(t) 



dxi 



j=i 



dxi 



vi=i 



pa(t) + pC.giz) J2 f i^r^' + 9{z) E 7./ (^)"^'^' / (^) ^. + P^.(0 



i=i 



i=i 



(39) 
(40) 

(41) 
(42) 



pCi\i + g{z) E / {zr-' + 9{z) E 7.7 i^r-' f {z) 



(43) 



-0, 



(44) 



where we require the following ordinary differential equations: 



a,{t) ^ F,{t) + C^i, fori = 1,2,. ..TV, 
a(0) = flo, a(0) = ai, 

e + 9{z) EU f i^r-' + 9{z) ET=i 7./ i^r^' f (^) = 0, 

5(2) > 0, for z G (—00,00). 



(45) 
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For 7i = 1, we have, 



p(??+|:-£)+|-K^+^^g^^M+^^^w (4^) 



^^('^ + ^ (^'"'^'^^) + ^ ^'"'^'^ ^ ^ ^"^'^ + ^^'(*) ^^^) 



^ -"^W + (^5 (^)) /(^) + ff(^)^/(^) (48) 

^ff (z) j C. ^ / (z)-^ + g{z)— Y. f i^r' + Pm) (49) 

: p-d[t) + pag{z)f (z) + g{z)f (z) C, (50) 

m in 

pC.g{z) ^ / {zY^'' + g{z) ^ 7,/ (^)''^'~' / (^) ^ + P^. W (51) 
i=2 3=2 

pcAi + pC,g{z) + 9{z)j^f^ + 5(^) E / (^)''~' + a{-A E ^J-/ (^)''"' / (^) [ (52) 

: 0, (53) 



with 



e + pC,g{z) + g{z)f§^C, + g{z) Y.^ f (^)'^"' + 9{z) Y.7=2 7./ (^)"'"' / (^) = 0, 

g{z) > 0, for z G (—00, cxd). 

The proof is completed. ■ 

With similar analysis, we can derive the corresponding theorem for the Euler-Poisson equations. 
Proof of Theorem [2l We can use the previous lemmas again to handle the mass equation and 
entropy equation again. Then, we can definite the potential function as: 

N 

$(i,f) = /(z) + ^d,(t)x„ (55) 



i=l 



N 

where z := T, Ci {xi — ai{t)) . 

i—l 

We differentiate ([55)) twice to obtain: 



A<J>(i,f) = /(z)ECf = a(iV)p, (56) 
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where in our solution 



N 



a{N) 



(57) 



For the i-th momentum equation ([51)2, we have: 



N 



dui s-^ duj \ d 

fe=i / ""'' \ j=i 

d 



d_ 
' dxi 



M^+E"^-a^ +^h^E^^M+^^7^'^+^^^w 



AT 



pau{t) + -^ e'"^(^) Y. h^V' + 5p-^ f{z) + 5]d.(i)x. + pF,{t) 



j=i 



ax,: 



pai 



771 t-j 777 

(t) + C,g (z) ^ / (z)^^- + g{z)— ^ / (z)^^ + 6p [j{z)C. + d.(t)) + pF.,{t) 



j=i 



j=i 



= pa.{t) + pGg(z) ^ / {z^ + 5(;2) E ^^-Z (^)'''"' / (^) ^^ 

+ (5p(/(z)G+d,(i))+pi^^(t) 



p^« u+ ^5wE/(^r"' + ^ffWE7,/(^r^' / (^) + ^^/w 



N 

7—1 



i=i 



7 = 1 



J = l 



0, 



(58) 

(59) 
(60) 
(61) 
(62) 

(63) 
(64) 



where we require the following ordinary differential equations: 



a,{t) = F,{t) + C,i + M,(t), for i = 1, 2, ...TV, 
a(0) = flo, d(0) = oi, 



e + ^¥^9{z) E™ 1 / (^)^^-^ + ^f^5(^) E^i 7,/ (^)"'^' / {z) + 5/(z) = 0, 



7.-2 



(65) 



z=l ^ 






(7(z) > 0, for z G (—00,00) 
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For 7i = 1, we get, 



fc=i / i=i 



d 



d 

dxi 



i=2 



+ ^P^ (fiz)+J2'^^it)x^ +pF,{t) 

= paS)+C,g{z)f{z)+g{z)-^J{z) 

+ Q<7(z)5]/(zp+5W— ^/(zp+5p(/(z)C. + d.(i))+pi^Ki) 

= pa,(i) + a.g(z)/(z) + 5(z)7 (z) C,+ 

■m m 

pa:5WE/(^)^' +-9WE^^/(^)^'"' 'fi^)C^ + Sp(^f{z)a+d,{t)) +pF,{t) 

j=2 3=2 



= pC^ < 



^+4^9i^) + 4^9i^)'i^ 



i=l ^ 






a(N) 



x7i-l , ajN] 



+^f^5(^) ET=2 f i^r + 4^ff(^) E" 2 7,/ (^r^^ / (^) + sfiz) 



i=l * 






(66) 
(67) 



(69) 



(70) 
(71) 

\ (72) 



0, 



(73) 



where we require the foUowing ordinary differential equations: 



e + 4^giz) + 4^9{z)il^ + 4^9iz) E"U f {^r~' 



s c^ 



i = l ' 



i = l ' 



a(Af) 



+^^ff(^) Er=2 7,7 (^)"^"^ / (^) + sm = 0, 



g(z) > 0, for z e (— oo, cx)). 



(74) 



The proof is completed. 



Remark 6 The ordinary differential equations (J 



^d^{t) = F,{t) + C^^ + rf,(<), /or « = 1, 2, ...A^, 
a(0) = ao, a(0) = ai, 



(75) 



are solved by 



at{t) 



JO 



(i^j(77) + di{ri)) drjds H ^- h ait + gq. 



(76) 
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Remark 7 In fact, our solutions, hl'2\) and \1^) . can he easily extended to the systems, with the 
generalized multiple and nonlinear, damping and Gamma pressure function: 

= -P^^^ + pF.it), for I = 1,2,. ..N, 
w/iere /3i,/32, ■•••,/3n; Ai, A2, .... ,Xrny and pi,p2, ...pn > 0, are constants. 
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